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A design study was conducted to investigate the effect shell buckling knockdown factor 
(SBKF), internal pressure and aluminum alloy material selection on the structural weight of 
stiffened cylindrical shells. Two structural optimization codes were used for the design study 
to determine the optimum minimum-weight design for a series of design cases, and included 
an in-house developed genetic algorithm (GA) code and PANDA2. Each design case specified 
a unique set of geometry, material, knockdown factor combinations and loads. The resulting 
designs were examined and compared to determine the effects of SBKF, internal pressure 
and material selection on the acreage design weight and controlling failure mode. This 
design study shows that use of less conservative SBKF values, including internal pressure, 
and proper selection of material alloy can result in significant weight savings for stiffened 
cylinders. In particular, buckling-critical cylinders with integrally machined stiffener 
construction can benefit from the use of thicker plate material that enables taller stiffeners, 
even when the stiffness, strength and density properties of these materials appear to be 
inferior. 
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I. Introduction 

M any components of launch vehicles consist of stiffened cylindrical structures, and these components may or 
may not be pressurized. Typical design practice for these cylinders uses conservative shell buckling 
knockdown factor (SBKF) values that are used to account for the differences that have been observed between test 
and theory. Recommended SBKF values for typical launch vehicle structure can be obtained from several sources, 
including references 1-3. Unfortunately, the pedigree of the test data, that typically dates from the 1930s to 1960s, is 
not well known, and in some cases (see ref. 3) does not account for pressure stabilization or length. As a result, 
conservative designs are developed that lead to weight penalty in a structure that is weight critical. In order to 
develop more efficient and effective launch vehicles, it is desirable to update the design process by developing less 
conservative SBKF values, and to account for pressure stabilization and for a variety of material and constraint 
types. 

As a first step in determining SBKF technology development priorities, a design study was conducted to 
investigate the portions of the typical design approach that can be improved to yield more structurally efficient 
cylinder designs for launch vehicle structure. The study used both an in-house developed genetic algorithm (GA) 
code and PANDA2 4 to examine unpressurized and pressurized cylinder designs. SBKF values were varied for both 
sets of designs, and pressure stabilization was also included. Additionally, inclusion of material constraints for 
related aluminum alloys demonstrates the importance of taking material size limitations into account in the design 
process for these cylinders. Optimum area weights for the study designs are examined and compared to show that 
use of less conservative SBKF values, including internal pressure, and proper selection of material alloy can result in 
significant weight savings for stiffened cylinders. 

II. Physical Description of the Structure 

Various cylinders with integrally machined orthogrid wall construction were considered in this study. Each 
cylinder section represents the uniform acreage design. No cutouts, weld land or additional structural details were 
accounted for in the design process. Definitions of the orthogrid geometry parameters are shown in Figure 1. The 
study considered three aluminum alloys that are defined in the following section. Designs were chosen that were 
considered to be typical of current launch vehicle structure. Unpressurized cylinders in this study had 200-inch 
radius and pressurized cylinders had 100-inch radius. Cylinder length (L) to radius ratios (L/R) of 0.5 and 0.75 for 
unpressurized cylinders and L/R of 1.5 for pressurized cylinders were used in this study. The cylinders were 
assumed to be simply supported on both ends. Applied loads are described subsequently, and consisted of uniform 
axial line load and internal pressure. 
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Figure 1 : Orthogrid geometry parameters. 


Material Properties 

The orthogrid cylinder designs utilized aluminum alloys whose properties were obtained from manufacturer data 
sheets 5 ' 6 or MMPDS-04 7 . Room temperature properties were used throughout the design study. Material information 
used in the design study for the three alloys used is shown in Table 1. All three alloys, 2195-T8R78 (2195), 2050- 
T84 (2050) and 2219-T851 (2219), were use in the unpressurized designs, while only alloy 2195 was used in the 
pressurized designs. 

Table 1 : Aluminum alloy material properties. 


Property 

2195-T8R78 3 

2050-T84 6 

2219-T851 7 

E (Msi) 

11.0 

11.1 

10.6 

a y (ksi) 

73.0 

71.0 

51.0 

a u (ksi) 

81.0 

75.0 

66.0 

p (lb/in J ) 

0.098 

0.098 

0.103 

V 

0.33 

0.33 

0.33 

Max. Plate Thickness (in.) 

1.8 

4.8 

5.8 


Loads 

Several load cases were considered for each of the stiffened cylinder designs, and include axial loads and/or 
pressure loads. The total axial load used for each load case represents the sum of an applied axial load and the 
maximum compressive contribution of an applied bending moment, and was assumed to be uniform around the 
cylinder circumference. This simplification was adopted as it is considered to represent the approach that is used in 
preliminary design. The total axial load is thereby defined by: 

N x = N xaxiaI + max comp^ N x<moment ) 

Load cases considered in this study include those that would represent proof pressure, pre-launch, flight and a 
fail-safe pressure-loss load case. Proof pressure is a load case whereby a pressurized tank component must sustain a 
prescribed pressure without yielding. The pre-launch load case represents the load while a launch vehicle is on the 
launch pad prior to applying the ullage pressure in tank components, and is subjected to wind loads and structural 
weight. The flight load case represents the worst loads experienced during the flight profile, and for unpressurized 
cylinders this was the only load case applied. Lastly, the pressure-loss load case represents a hypothetical fail-safe 
flight load case where the ullage pressure is reduced to some prescribed level and global buckling is not permitted 
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prior to limit load (SF=1.1). This approach is used in designs that consider pressure stabilization. The proof load 
case was used prior to optimization to determine the minimum skin thickness for pressurized cylinders by requiring 
that the skin sustain the proof hoop stress without yielding. Mathematically, this is given by: 

_ P proof R 

1 skin, min 


where t skin ^ n is the minimum skin thickness, P pr0 of is the proof pressure, R is the cylinder radius and a v is the yield 
stress of the material. In this study, P pmo f = 66 psi, resulting in a minimum skin thickness of 0.09 inches 

Unpressurized cylinders were optimized under flight loads of N x = -1800 and -2000 lb/in. Load case definitions 
for the pressurized cylinders are shown in Table 2. In the pressure-loss case Pi oss = 0.0 psi, optimization was 
performed for three load set combinations; 1) P rdief = P hoop-stabmty = 0.0 psi, 2) P axmMi( , f = 20.0 psi and Ph oop - s tabiiity = 
0.0 psi, and 3) P m mi-rdwf = P 'hoop-stability = 20.0 psi. Maximum and ullage pressures were assumed to be P max = 50.0 psi 
and Puiiage = 30.0 psi in the flight load case for stress calculations in the pressurized cylinders. Only aluminum alloy 
2195-T8R78 was used in the pressurized cylinder portion of the study. 

Table 2: Load case definitions for pressurized cylinders. 


Load Case 

Nx (lb/in) 

P axial-relief (psi) 

Phoop-stabilitv (psi) 

Pioss (psi) 

Prelaunch 

-800 

0.0 

0.0 

N/A 

Flight 

-3200 

0.0, 20.0 

0.0, 20.0 

N/A 

Pressure-loss 

-3200 

p 

1 loss 

p 

1 loss 

0.0, 1.0, 2.0, 4.0, 
6.0, 10.0, 20.0 


III. Design Requirements and Criteria 


Structural Design 

Designs were developed using two different optimization codes to determine the sizing of the various orthogrid 
geometry parameters shown in Figure 1. Numerous design constraints were assigned that are categorized as failure 
constraints and geometry constraints. Failure constraints included; 

• Global buckling 

• Local buckling (see ref. 8) 

• Limit strength (unidirectional in axial and circumferential directions) 

• Limit stress (bi-axial) 

• Ultimate strength (unidirectional in axial and circumferential directions) 

• Ultimate stress (bi-axial) 

• Stiffener buckling (PANDA2 only) 

• Stiffener crippling (see ref. 9, GA only) 

• Crippling/buckling interaction (see ref. 9, GA only) 

Geometry constraints included; 

• Shell thickness (in.) (t skin>min < t s < 0.40 in.) 

o Pressurized: t skinjni „ derived from proof pressure, 0.09 in. 

o Unpressurized: t skin min derived from minimum manufacturing gauge of 0.065 in. 

• Stringer spacing (2.0 in. < d str < 16.0 in.) 

• Stringer thickness (in.) (0.065 in. < t str <1.0 in.) 

• Ring spacing (2.0 in. < d str < 50.0 in.) 

• Ring thickness (in.) (0.065 in. < t ring < 1.0 in.) 

• Ring aspect ratio (0.0 < AR lirlg < 12.0) 

o Manufacturing constraint due to forming process 

• Stringer height = Ring height 

• Minimum stringers per circumferential wave number, n=10 

o GA only, to justify the use of smeared stiffener approximation 

• Maximum plate stock thickness (material alloy dependent, see Table 1) 

4 

American Institute of Aeronautics and Astronautics 


Tracking #210774 


Only internal stiffeners were considered in the current study, and the height of the stiffeners is defined as the 
distance from the outer mold line to the top of the stiffener. The minimum and maximum limits for the stringer and 
ring spacing were arbitrarily chosen for the design study. The ring aspect ratio, AR ring , is the ring height divided by 
the ring thickness, h ring /t ring ,, with the maximum limit value of 12.0 being set by manufacturing constraints. 


Safety and Knockdown Factors 

Relevant safety factors and knockdown factors for this type of structure were used in the design process. The 
safety factors apply to all design loads and include limit stress, ultimate stress, global buckling and local buckling. 
Knockdown factors for global buckling are applied to account for reduction in predicted buckling load due to 
geometric imperfections. The baseline value for the shell buckling knockdown factor for stiffened shells, y = 0.65, 
was obtained from the Isogrid Design Handbook 3 . Larger magnitude shell buckling knockdown factors are used in 
pressurized cylinders to account for the fact that internal pressure reduces the sensitivity of the buckling load to 
initial geometric imperfections. This increase, Ay, is dependent upon the internal pressure of the cylinder. The Ay 
values used for the pressurized cylinders were developed from an imperfection sensitivity study for pressurized 
cylinders using the STAGS 10 finite element code. These analysis-based Ay values are smaller in magnitude than 
would be obtained using the method described in SP-8007 1 , however, they represent a conservative estimate of the 
potential benefit of reduced imperfection sensitivity associated with pressurized cylinders. The study was begun 
using the baseline shell buckling knockdown factor of y — 0.65. This value was then increased to y = 0.75 and y = 
0.85. The SBKF for pressurized load cases was then increased by Ay = 0.14 when pressure stabilization was 
included for the cylinders being studied. The SBKF for the pressure-loss load case was linearly interpolated between 
the unpressurized SBKF and the pressure-stabilized SBKF. Mathematically, this is represented by: 


Y loss 


y + Ay 


p 

y hoop- stability 


\ 

/ 


A local buckling knockdown factor of 0.9 was chosen, and is considered to be very conservative. 

The safety factors applied to the design cases in this study are as follows; 

• Limit =1.1 

• Ultimate =1.4 

• Global buckling =1.4 

o Reduced to 1.1 or 1 .2 in pressure-loss case for pressurized cylinders 

• Local buckling =1.25 

Margin Calculations 

Margins of safety (MS) were calculated based upon allowable values or design requirements. A generic form of 
the MS calculation equation is: 


. Specified Value 

MS = — - — = 1 

Calculated Value 

For a stress margin of safety, the calculation was performed using the following equation: 


MS = 


° allow 

SF(o calc ) 


-1 


where SF represents the applicable safety factor, a a n ow is the allowable stress (a y for limit load and o u for ultimate 
load), and a ca i c is the calculated stress. The buckling margin of safety calculations include the appropriate 
knockdown factor as given in the following equation: 


MS = 


y(AC fc ) 

sf(n sp<x ) 
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where y is the knockdown factor, N c c r aIc is the predicted critical buckling line load for a given design, and N spec is the 
specified design line load. The design study required that all calculated MS be greater than or equal to zero. 

IV. Design Approach 

Designs for the cylinders were obtained by using two different optimization design codes, a genetic algorithm 
(GA) based code and PANDA2 4 . The GA code was developed in-house by linking a GA with BOSOR4 11 , where 
BOSOR4 was used to determine global buckling loads and mode shapes. Stiffeners, rings and stringers for the 
designs were approximated using a smeared stiffener approach in the global buckling analyses. Additional constraint 
calculations in the GA code, such as allowable stress, local buckling and crippling, are similar to those used in other 
commercially-available design codes and handbooks. For the GA, ten random number seeds were used to initiate ten 
optimization analyses for each design, and the best solution was chosen as the minimum weight design that satisfied 
all the design constraints. PANDA2 is a panel optimization code that also uses BOSOR4. For this study, a single 
PANDA2 optimization analysis was used, with three or four superopt 4 executions being conducted for each design. 
It should be noted that PANDA2 has built-in constraints additional to those discussed in Section III. 

PANDA2 permits only one material stress allowable, so an additional load case was required to examine both 
yield and ultimate stress constraints. The additional PANDA2 stress load case was only applied to the flight load 
case since it was assumed that the other load cases resulted in equal or lower stresses, and therefore, the role of 
stress in the design was accounted for in the flight load case. At the same time, the buckling knockdown factors for 
this additional load case were increased because buckling was already being addressed by the original flight load 
case. The separate yield and ultimate stress analysis capability is included for each load case in the GA code by 
accepting both yield and ultimate stress allowable values. 

V. Results 

Area weight results for unpressurized orthogrid cylinders are shown in Figures 2 - 5. It is seen that the results for 
the GA and PANDA2 2195 designs are in very good agreement for area weight and weight reduction trend. In 
contrast, the designs for 2050 and 2219 show good agreement for weight reduction trend, but show significant 
difference in area weights. This difference likely arises from the fact that PANDA2 has additional constraints and 
additional conservatism applied to the design compared to the GA code. For example, for the 2219 design weights 
shown in Figure 2 with y = 0.65, N x = -1800 lb/in., and L/R = 0.5, the GA and PANDA2 designs have area weights 
of 2.05 and 2.68 lb/ft 2 , respectively. Examination of the constraints show that both consider global buckling and 
local buckling as design drivers, with both having MS of near zero. The PANDA2 design has an additional 
constraint, lateral buckling on the stiffener, that also has an MS of nearly zero. This constraint is not included in the 
GA code. As a result, the stiffener height in the GA design is 2.23 in. as compared to 1.87 in. for the PANDA2 
design. Similarly, the stringer aspect ratio is much larger for the GA design than the PANDA2 design. Therefore, the 
GA code allows for a stiffener that is much more efficient when it comes to global buckling, but violates a constraint 
that is included in PANDA2. This discovery suggests that these constraints should be investigated further and the 
GA code modified as necessary. However, the current study is concerned with the trends in weight reduction, and 
the designs from the two codes show very good agreement in this respect. 

It is seen in the Figures 2-5 that for the 2195 designs, there is significant weight reduction when the SBKF is 
increased from 0.65 to 0.85, with area weight reductions on the order of 18 to 24 percent. However, compared to the 
2195 designs, area weight reductions for the 2050 and 2219 alloys are significantly smaller, on the order of 4 to 13 
percent for the same range of SBKF values, and are much more uniform across the range. These trends are due to 
the fact that the smaller plate stock thickness of the 2195 alloy imparts a weight penalty resulting from the material- 
induced stiffener height constraint. As a result, the designs use the entire available plate thickness, and the 
thicknesses of the stringers and rings are required to be much thicker in order to provide the required bending 
stiffness necessary to resist global buckling. Contrary to this, the designs for the 2050 and 2219 alloys take 
advantage of the increased material thickness by developing taller stiffeners, but never approach requiring the entire 
plate thickness. 

The GA results in Figures 2-5 show clear separation in the area weight of the designs for the three alloys. 
However, the PANDA2 results suggest that the 2050 designs result in minimum weight for all SBKF values in this 
study, but for 2195 and 2219 the minimum weight cylinder is dependent upon the SBKF value. Consider the 2219 
design weights shown in Figure 3 with y = 0.65, N x = -2000 lb/in., and L/R = 0.5, where the GA and PANDA2 
designs have area weights of 2.15 and 2.86 lb/ft 2 , respectively, while the corresponding 2195 designs have GA and 
PANDA2 area weights of 2.87 and 3.00 lb/ft 2 , respectively. Despite 2219 having a density that is 5.1 percent higher 
than 2195, while at the same time having a stiffness (modulus) and a strength allowable (yield) that are 3.6 and 30 
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percent lower than 2195, respectively, both codes yield 2219 designs that result in reduced weight compared to 
2195. Weight savings predicted by the GA and PANDA2 for these designs are 25 percent and 4.7 percent, 
respectively. The reduced savings for the PANDA2 designs is attributed to the aforementioned additional constraints 
included in PANDA2. Next, for N x = -2000 lb/in. and L/R = 0.5, but with y increased from 0.65 to 0.85, PANDA2 
yields a 2195 design that is 7.2 percent lower than the corresponding 2219 design. This switch in optimum material 
results from the lower stiffness of the 2219 having a greater impact on the stiffener lateral buckling constraint in 
PANDA2. Therefore, depending on the design code used and the active constraints for the design, it is possible that 
material selection should become a design variable in order to obtain minimum weight designs. 
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Figure 3: Area weight of unpressurized, orthogrid cylinder; N x = -2000 lb/in., L/R = 0.5. 
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Figure 4: Area weight of unpressurized, orthogrid cylinder; N x = -1800 lb/in., L/R = 0.75. 
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Figure 5: Area weight of unpressurized, orthogrid cylinder; N x = -2000 lb/in., L/R = 0.75. 

Area weight results for pressurized, orthogrid cylinders fabricated from 2195 alloy are shown in Figures 6-8. 
Results for the pressure-loss load case with global buckling SF = 1.1 and 1.2 are shown in Figures 6 and 7, 
respectively. These figures also show the area weight for when no pressure relief or stabilization is included in the 
design for the typical approach where y = 0.65 (black circles). Also shown is the case where pressure relief is 
included but pressure stabilization and pressure-loss load case are not, and where y = 0.65 (black squares). When 
pressure stabilization is included, the pressure-loss load case becomes the design driver up until the reduced internal 
pressure reaches a certain value, as indicated by an abrupt change in slope on the plots. As seen in the figures, the 
pressure value at which the controlling load case transitions from the pressure-loss load case to the flight load case is 
a function of global buckling SF in the pressure-loss load case. Above this transition pressure, the pressure-loss load 
case no longer affects the design. For example, changing the global buckling SF value from 1.1 to 1.2 changes this 
transition pressure value from about 12.5 psi to about 14 psi. 

Maximum weight savings can be obtained by including pressure relief and stabilization, and by omitting the fail- 
safe pressure-loss load case, with area weight reductions of 36 to 41 percent for SBKF values from 0.65 to 0.85. 
When including the pressure-loss load case, where the internal pressure is reduced to 0.0 psi, the weight reductions 
range from 8.5 to 24 percent. Next, a comparison of the area weights as a function of the global buckling SF in the 
pressure-loss load case is shown in Figure 8 for y = 0.65. It is clear that increasing the SF value will result in a 
reduction in weight savings. Judicious choice of global buckling SF and reduced internal pressure in a pressure-loss 
load case can result in significant weight reduction when compared to an unpressurized design, while providing for 
various fail safe options. 
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Figure 6: Effect of pressure-loss load case P hss value on area weight of pressurized cylinder, PANDA2 
results; pressure-loss load case global buckling SF= 1.1. 
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Figure 7: Effect of pressure-loss load case Pi oss value on area weight of pressurized cylinder, PANDA2 
results; pressure-loss load caseglobal buckling SF= 1.2. 
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Figure 8: Effect of pressure-loss load case global buckling 5F value on area weight of pressurized 

cylinder, PANDA2 results; y = 0.65. 

VI. Conclusion 

The stiffened cylinder designs in the current study were examined to determine the effects of the SBKF, internal 
pressure and material selection on the optimized design weight. Weight comparison was based upon a traditional 
approach without pressure relief or stabilization. In particular, it was shown that SBKF and internal pressure can 
have a significant effect on the design weight of buckling-critical designs. For example, weight reductions 
approaching 41 percent were achieved by using a higher SBKF, and by including hoop stabilization due to internal 
pressure, which accounts for the stabilizing hoop stress and an additional increase in SBKF Including a pressure-loss 
load case to provide a fail-safe condition results in a lower weight savings. This weight savings is dependent upon 
the global buckling SF and reduced internal pressure used in the pressure-loss load case. Flowever, based on the 
typical design approach, including a pressure-loss load case can still yield significant weight savings. Additionally, 
for unpressurized cylinders, selection of a material that is manufactured in thicker stock can provide a significant 
weight savings. This weight savings is expected to be applicable to pressurized cylinders, as well. These 
observations lead to three general conclusions: 

1. Reasonable increases in SBKF, from y = 0.65 to y = 0.75 or 0.85, can result in lower cylinder design 
weights. Therefore, it would be advantageous to develop less conservative SBKF values rather than 
relying on the conservative values that are typically used. 

2. Allowing for hoop stabilization can result in lower cylinder design weights. Flowever, if additional safety 
is required in a design, a fail-safe pressure-loss load case can be employed and still yield significant 
weight savings. 

3. It is important to take into account the material dimensions in the design of cylinders. Due to 
manufacturing processes, various alloys may be available in widely different plate stock thicknesses. 
Therefore, it is important to understand the design drivers. For example, when the design is buckling 
critical, use of a material that comes in thicker stock can result in significant weight savings. This weight 
savings arises from the potential for increased stiffener height, which yields more effective stiffeners. The 
potential weight savings can even be achieved when the thicker plate stock material has what could be 
deemed as less desirable properties, such as increased density and lower stiffness and strength. 
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